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We evaluate the accuracy of traditional estimators often used to extract TV — » A quadrupole 
ratios from cross section angular distributions for pion electroproduction. We find that neither 
Mi+ dominance nor I < 1 truncation is sufficiently accurate for this purpose. Truncation errors 
are especially large for Rem, for which it is also essential to perform Rosenbluth separation. The 
accuracy of similar truncated Legendre analyses for Eq+, So+, and especially Mi_ is even worse. 
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Historically the most important indications of deformation of low-lying baryons have been the quadrupole ratios for 
electromagnetic excitation of the N — > A(1232) transition. Magnetic dipole excitation dominates and is represented 
by the M\ + multipole amplitude while nonzero values for the electric and scalar (longitudinal) multipoles, E\+ and 
Si+, arise either from nonsphcrical contributions to the wave functions or from higher-order dynamical contributions 
to the electromagnetic transition. The quadrupole ratios are defined as 

Rem = Re^i (la) 

Mi + 

Rsm = Re^±- (lb) 

M 1+ 

evaluated at the physical mass of the resonance, W = Ma ~ 1.232 GeV. Determination of quadrupole ratios for 
isospin-3/2 amplitudes requires measurements of two charge states, such as pir° and nir + . Complex multipole am- 
plitudes have been deduced for Q 2 — using polarization data for pion photoproduction but few experiments 
for Q 2 > have provided sufficient information to perform an actual multipole analysis. Instead, most experimental 
determinations of N — > A transition form factors [2J, |3J, |J, |5j rely upon estimators derived from multipole expansions 
for the angular dependence of unpolarized cross sections using two simplifying assumptions: 1) only multipoles with 
£ < 1 contribute, which is described as sp truncation; and 2) only terms involving M\+ are retained, which is de- 
scribed as M 1+ dominance. Although the reliability of these assumptions has been questioned before, the improved 
kinematic completeness and statistical precision of modern experiments warrants re-examination of their accuracy. In 
this Brief Report, we consider the accuracy of traditional quadrupole estimators for Q 2 ~ 1 (GeV/c) 2 where a nearly 
model-independent multipole analysis of recoil-polarization response functions for the p(e, e'p)ir° reaction disagrees 
appreciably with the traditional Legendre analysis of the cross section data [f| Q • 

The unpolarized differential cross section for j v N — > Nn in the irN center of momentum frame takes the form 

da 



vo [csRl + Rt + \/2es(l + e)R LT sin 9 cos (j) + eR T r sin 2 9 cos 2 cf>) (2) 



where vq is a phase-space factor, e is the transverse polarization of the virtual photon, £5 = eQ 2 /q 2 , and (9,<fi) are 
polar and azimuthal pion angles relative to the q vector and the electron scattering plane. The response functions 
can be expanded in Legendre series 

00 

R X = J2 A n P n( COS0 ) ( 3 ) 
n=0 

where A € {L,T, LT,TT}. The expansion coefficients, A„, are functions of (W,Q 2 ) that can be fit to the angular 
distribution of the differential cross section. Each of those coefficients can in turn be expressed as a multipole expansion 
containing terms of the form Rei?^±C| ± where B,C £ {M, E, S} are magnetic, electric, or scalar multipole amplitudes 
for specified I and j = £ ± 1/2. In principle, experimentally determined Legendre coefficients include contributions 
from arbitrarily large £ and are not limited either by sp truncation or by M14- dominance. 

Truncated multipole expansions of the Legendre coefficients used in quadrupole estimators are given in Eq. 
where the contributions that satisfy M\+ dominance are listed first and where the remaining terms include some of the 
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TABLE I: Complexity of multipole expansions of Legendre coefficients. For each Legendre coefficient used for quadrupole 
estimators we show the number of independent terms of the form Re(ab*), where a and b are multipole amplitudes with 
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lowest multipolarity contributions of other types but are not necessarily arranged in order of numerical importance. 

Aq = |S 0+ | 2 + 8|S 1+ | 2 + |S 1 _| 2 + 8|S 2 _| 2 + 27|S 2+ | 2 + ... (4a) 
Al = 2|M 1+ | 2 + |So + | 2 + |A/ 1 _| 2 + 6|^i + | 2 ]+6|Af 2 -| 2 + 2|S 2 _| 2 + 9|il/ 2+ | 2 + 18|S 2+ | 2 + ... (4b) 

Al T = ~|Af 1+ | 2 - Re[M* + (3£i+ + 3M X _ + 12M 3 _ + 3£ 3 _ + 2M 3+ + 10£ 3+ )] (4c) 
+ ^-\E 1+ \ 2 + ||£7 2 _| 2 + 24|£ 2+ | 2 - ||A/ 2 _| 2 - 12|M 2+ | 2 

+ Re[-3^ + (^ 2 - + M 2 _ - M 2+ + E 2+ ) + E* + (3M^ - 21£ 3 _ - 12M 3 _ + 12M 3+ ) 
+ M*_(3E 3 _ + 3Af 3 _ + WE 3+ - 10M 3+ ) + ...] 
A\ T = 3Re[M 1 * + (25i+ - 3S 3 _ + 45 3+ ) + S£ + (E 2 - - M 2 _ + Af 2+ - AE 2+ ) (4d) 
- S *+(2^2- - 3S 2+ ) - 2E* 1+ (S 1+ + Sj.-) - 2M 1 _5 1 % + ...] 
216 

^2 = 8|5 1+ | 2 + 8|5 2 _| 2 + — |5 2+ | 2 +Re[5 * + (85 2 _ + 185 2+ ) + 85 1+ 5 1 *_ + ...] (4e) 

Al = -|M 1+ | 2 +Rc[Afr + (6S 1+ -2Af 1 _ + yAf 3 _+6^ 3 _ + ^M 3+ )] (4f) 

+ 3|£ 1+ | 2 - |£ 2 _| 2 + ^|£ 2+ | 2 + 3|Af 2 _| 2 + y |A/ 2+ | 2 

+ Re[.E * + (2£ 2 _ - 6M 2 _ + 6M 2+ + f 2£ 2+ ) - 6Mi_-E? + + ...] 

Table [I] shows that the number of independent terms in the multipole expansions of these Legendre coefficients 
increases very rapidly with the maximum £ permitted. Complete expansions for i" max < 6 can be found in Ref. Q but 
as £ max increases they quickly become too unwieldy to display here or to use in practical applications. The Legendre 
coefficients are usually obtained by numerical integration of response functions against Legendre functions instead of 
by these algebraic formulas, but both methods do agree. 

The assumption of AT 1+ dominance omits any terms that do not involve Af 1+ , which strongly inhibits the prolifer- 
ation of terms but is not sufficient in itself to extract quadrupole ratios from cross section data. Combined with sp 
truncation, these expansions reduce to 

Aq « (5a) 
Al a 2|Af 1+ | 2 (5b) 

Al T a -Re[(|Mi++3£7i++3Mi_)Af 1 * + ] (5c) 

A[ T a 6Re[Si+Af* + ] (5d) 

A 2 « (5e) 
Al a Re[(-Af 1+ + 6^1+ - 2Af 1 _)M 1 * + ] (5f) 

Thus, one obtains the traditional quadrupole estimators 



3( Al + eAh)-2A l T 
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A LT 

RsM = mTTAf) (6b) 
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TABLE II: Convergence of multipole expansions of Legendre coefficients and quadrupole estimators. Multipole amplitudes 
from MAID2003 for the pn channel were used for W = 1.232 GeV and Q 2 = 1.0 (GeV/c) 2 . Legendre coefficients are in units 
of^fc) 1 / 2 . 

















e 


= 0.95 


c 


= 


p 

max 




A T 




aTT 
Aq 


Ak 


Al 


Jem 


fsM 


I EM 


fsM 


l 


0.3339 


7.599 


-1.483 


-5.220 


0.141 


-3.769 


0.300 


0.938 


0.581 


0.977 


2 


0.3377 


7.624 


-1.395 


-5.156 


0.102 


-3.916 


0.736 


0.879 


0.960 


0.916 


3 


0.3384 


7.628 


-1.323 


-5.137 


0.079 


-3.871 


0.714 


0.833 


0.894 


0.868 


4 


0.3384 


7.628 


-1.315 


-5.130 


0.081 


-3.859 


0.696 


0.828 


0.878 


0.863 


5 


0.3384 


7.628 


-1.308 


-5.130 


0.081 


-3.857 


0.693 


0.824 


0.876 


0.859 



where is included because most experiments have not used Rosenbluth separation to isolate A^. When Rosenbluth 
separation is available, one can simply use e — > in Eq. J()J|. Therefore, it is useful to define the accuracy parameters 

f (t A 1 f ^ + eA!;)-2Ar \ „ Rem , . 

/sM(Cax ' £) = flb (w^L,_ " iff (7b) 

where asymptotic equality refers to the limit f max — > oo. Despite their appealing simplicity, it is clear that many 
contributions are omitted and the accuracy of the traditional estimators is a numerical issue that can be addressed 
either theoretically using model calculations or experimentally using additional polarization measurements to extract 
complex multipole amplitudes directly. 

The convergence of these expansions is evaluated in Table ITU using pn multipole amplitudes for W — 1.232 GeV 
and Q 2 = 1.0 (GcV/c) 2 from MAID2003 HG3- First ' we observe that A„ contributions are not negligible: the 
contribution of Aq to the denominators of Eq. @) reduces the estimated quadrupole ratios by about 4% without 
Rosenbluth separation when e — ► 1. (Note that e = 0.949 at W = 1.232 GeV in Ref. @.) Even though A^ is much 
smaller, its effect upon Jem is even larger because the strong cancellation between Aq T and Aj + tA\ amplifies 
the dependence on e. Therefore, the assumption of Mi+ dominance is not sufficiently accurate to measure Rem 
without Rosenbluth separation. Even with Rosenbluth separation, one should not expect better than about 15% 
accuracy for either quadrupole ratio using the traditional Legendre analysis (see the bottom of last two columns of 
Table ILT|) . Second, it is clear that sp truncation is not valid either because contributions with I > 1 are not negligible. 
Cancellation between contributions to the numerator of Rem also amplifies truncation errors and convergence is 
not necessarily monotonic as £ max increases. Contributions to Aq and Aq are nonnegative, but the signs for other 
Legendre coefficients are mixed. While the magnitudes of multipole amplitudes for £ > 1 do tend to decrease, their 
coefficients in Eq. (@J tend to increase with £. Thus, convergence becomes a delicate numerical issue. 

Under the present conditions, we find that Re(Mi_£' 1 ' + ) is the most important contribution to Rem neglected by 
Mi_l_ dominance and is approximately —40% of the leading term. Thus, Mi+ dominance is not very accurate. The 
fact that Jem approaches 0.88 for e = is actually nothing more than a lucky conspiracy among the magnitudes and 
signs for a very large number of smaller terms, many of which are not especially small individually. However, most 
experiments omit Rosenbluth separation. Similarly, the second most important contribution to Rsm is Re (So +£"2-) 
but is only about 6% of the leading term; hence, /sm converges more rapidly. The details of this analysis are obviously 
model dependent, but qualitatively similar results are obtained for other models as well. Although Jem for e = 
is slightly closer to unity than /sm for the present analysis, the greater susceptibility of Rem to truncation errors 
through its reliance upon delicate cancellations suggests that the traditional Legendre analysis is intrinsically less 
reliable for Rem than for Rsm, with or without Rosenbluth separation. 

It is often argued that the traditional Legendre analysis should be more accurate for the isospin-3/2 channel than 
for the pir° reaction because the resonant multipoles should share a common phase and become pure imaginary at 
the physical mass, thereby suppressing background contributions. Leaving aside the propagation of errors involved in 
extracting isospin-3/2 amplitudes by combining two independent experiments, we can address the intrinsic accuracy 
of this analysis method using model calculations also. The convergence of the accuracy parameters for isospin- 
3/2 quadrupole ratios is examined in Table ITTT1 Again we find that Rosenbluth separation is required for Rem- 
Interestingly, /em deteriorates as f max increases and the final accuracy of Rem is worse for isospin-3/2 than for 
pir° even with e = 0. The cancellations are severe, the method is unstable, and calculations for Rem are highly 
model-dependent. 
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TABLE III: Convergence of quadmpole estimators for isospin-3/2. Multipole amplitudes from MAID2003 were used for 
W = 1.232 GeV and Q 2 = 1.0 (GeV/c) 2 . Legendre coefficients are in units of {pb) 1/2 . 
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FIG. 1: W dependence of quadrupole estimators for piv at Q 2 — 1.0 (GeV/c) 2 using MAID2003 multipoles. Solid curves show 
Rem and Rsm while dashed and dash-dotted curves show Rem and Rsm for I < 5 using e = and e = 0.9, respectively. The 
vertical dashed line denotes W = 1.232 GeV. 



Figure n compares traditional quadrupole estimators with Rem an d Rsm f° r MAID2003 at Q 2 = 1.0 (GeV/c) 2 . 
Ideally the estimators would be most accurate in the immediate vicinity of the physical mass, W = Ma ~ 1.232 GeV, 
but neither actually has that property. Rosenbluth separation is most important for Rem, but even with separation 
the residual error at Ma is significant at the level of experimental precision that is now possible. 

Similarly, Fig. [21 shows the Q 2 dependence for the accuracies of the traditional quadrupole estimators using 
MAID2003 p7r° multipole amplitudes for W — 1.232 GeV. Solid curves use e = 0, corresponding to Rosenbluth 
separation, while dashed curves use e = 0.9, typical of many experiments. Even with Rosenbluth separation, neither 
quadrupole estimator can be trusted to better than about 20% and their accuracy deteriorates at larger Q 2 as Mi + 
dominance breaks down. Therefore, truncation errors can seriously affect the Q 2 dependence of quadrupole amplitudes 
deduced from Legendre coefficients. The most extensive recent study of pit quadrupole ratios for 0.4 < Q 2 < 1.8 
(GeV/c) 2 used Eq. © without Rosenbluth separation 0. We do not advocate adjustment of such results using Fig. 
13 at least at this time, because the shapes of /em and Jsm are model dependent and MAID2003 does not describe 
all of the low- lying multipole amplitudes at Q 2 — 1 (GeV/c) 2 from Ref. sufficiently well to be confident of its 
predictions for the Q 2 dependencies of these ratios. Instead, we claim that accurate measurements of the quadrupole 
ratios require multipole analysis of both polarization and cross section data. 

Finally, other simple estimators 

RcE Q+ M* + w Al/2 (8a) 
ReSo+Mj^ « A^ T (8b) 
ReMi-M^ « -{2Al + 2Al T + Al)/% (8c) 

based upon A/i + dominance and sp truncation are sometimes quoted [21 1 1 1 1] - Note that Rosenbluth separation is 
required. However, using MAID2003 pir° multipoles at (W, Q 2 ) = (1.232, 1.0) with £ < 5, the ratios between the right- 
and left-hand sides of Eq. © are 1.74, —0.77, and 9.75. Most notably, the numerical contribution of ReA/i-M-j^ is 
only the fifth largest term in the multipole expansion of the specified combination of Legendre coefficients. Therefore, 
these estimators are worthless under these conditions. 

In summary, we have performed a detailed numerical analysis of truncation errors in quadrupole ratios deduced 
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FIG. 2: Accuracy of traditional quadrupole estimators for p7r° at W = 1.232 GeV using MAID2003 multipoles. Solid curves 
use e = and dashed curves use e = 0.9. 



from Legendre coefficients fit to cross section angular distributions. We find that neither M\ + dominance nor sp 
truncation is reliable and that one cannot expect better than 20% accuracy from this method. Truncation errors are 
especially important for Rem- Furthermore, accurate results for Rem also require Rosenbluth separation, which was 
not performed in recent studies of the Q 2 dependence of the quadrupole ratios. The accuracy of truncated Legendre 
analyses of Eq+, Sq+, and especially Mi_ is even worse. Polarization data for pion electroproduction are needed 
to perform nearly model-independent multipole analyses that provide complex amplitudes and do not depend upon 
unjustifiable truncation schemes. 
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